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We investigate the one-dimensional propagation of waves in the Anderson localization regime, 
for a single-mode, surface disordered waveguide. We make use of both an analytical formulation 
and rigorous numerical simulation calculations. The occurrence of anomalously large transmission 
coefficients for given realizations and/or frequencies is studied, revealing huge field intensity concen- 
tration inside the disordered waveguide. The analytically predicted s-like dependence of the average 
intensity, being in good agreement with the numerical results for moderately long systems, fails to 
explain the intensity distribution observed deep in the localized regime. The average contribution 
to the field intensity from the resonances that are above a threshold transmission coefficient T c is a 
broad distribution with a large maximum at/near mid- waveguide, depending universally (for given 
T c ) on the ratio of the length of the disorder segment to the localization length, L/£. The same 
universality is observed in the spatial distribution of the intensity inside typical (non-resonant with 
respect to the transmission coefficient) realizations, presenting a s-like shape similar to that of the 
total average intensity for T c close to 1, which decays faster the lower is T c . Evidence is given of 
the self-averaging nature of the random quantity log[/(a;)]/a; ~ — 1/£. Higher-order moments of the 
intensity are also shown. 



PACS numbers: 72.15.Rn, 42.25.Dd, 72.10.Fk 

I. INTRODUCTION 

There are two well-known manifestations of strong lo- 
calization of classical waves in one-dimensional (ID) open 
disordered systems: exponentially small (with respect to 
the length of the system) transmission through typical 
(most probable) random realizations, and high trans- 
parency at rare (exponentially low-probable) ones. The 
high transparency is due to the so called stochastic reso- 
nances that are accompanied by large concentration (lo- 
calization) of energy in relatively small areas inside the 
system. It was shown in Ref. yjthat in a semi-infinite ran- 
dom medium the wave amplitude at the resonances can 
exceed (with nonzero probability) any given value. In the 
80s this phenomenon had been studied intensively as ap- 
plied to electrons, light, elastic, and acoustical wavesS* 3 - 
(see also Refs. IjlalaQand references therein). In the last 
few years, after a long hiatus, interest in stochastic reso- 
nances in random media has rekindled in the context of 
random lasing 8 910-11 , wherein resonances might play the 
role of effective confining cavities inducing lasing action 
when gain is introduced. 

We investigate the one-dimensional propagation of 
electromagnetic (EM) waves in the strong localiza- 
tion regime. In particular, the occurrence of anoma- 
lously large transmission coefficients for given realiza- 
tions and/or frequencies (resonant or quasi-transparent 
realizations) is studied, with emphasis on the field inten- 
sity distributions along the direction of propagation. For 
that purpose, we make use of both an analytical formu- 



lation and rigorous numerical simulation calculations. 



We consider a single-mode waveguide with randomly 
rough walls. This structure, being a typical example, of 
a one-dimensional disordered system, has the advantage 
that it can be easily prepared using standard equipment 
(microwave waveguides or fiber optics), and enables (un- 
like a random stack of dielectric layers) to directly mea- 
sure the wave field inside the structure. Similar multi- 
mode systems have been studied in recent years to in- 
vestigate various localization and transport phenomena 
appearing in the propagation of waves through disordered 
media 12 i 13 i 14 i 15 i 16 i 17 i 18 . 



Our numerical calculations exploit the invariant em- 
bedding equation formulation for a multi-mode surface- 
disordered waveguid e 12 ' 14 ! 19 , which we have extended to 
account for the field inside the disordered region. The 
numerical results are compared with analytical formulas 
obtained by using the invariant embedding method and 
averaging over rapid phase variations 4,19 . Both methods 
are described in Sec. [D] Local and average field inten- 
sities are presented Sees. 11111 and IIVI respectively; the 
conclusions drawn from them are summarized in Sec.lVl 
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II. SCATTERING MODEL 

A. Field distribution outside the disordered region: 
Reflection and transmission amplitudes 

The scattering geometry is depicted in Fig. We seek 
for solutions to the scalar Helmholtz equation in the form 
(outside the region < x < L): 

(x, r) = k m 1/2 Xm (r) e- lkmX t mn , x < 0, (la) 

m 

* n (x,r) = #° (x,r) 

+ fc m 1/2 Xm (r) e tkmX r mn , x > L, (lb) 



with 



*° (x,v) = k-V 2 Xn (r)e- ik ^. 



(lc) 



The indexes "m, n" correspond to the outgoing and in- 
coming modes, respectively. \ n (r) are the eigenfunctions 
of the transverse wave equation, characterized by trans- 
verse momentum K n , so that the longitudinal wavevector 
component k n (along the propagation direction) is 



kn = [(^/c) 2 - kI 



21 1/2 



(2) 



with uj being the wave frequency. 

We consider the Dirichlet boundary condition on a 
slightly perturbed waveguide surface, C denoting the ran- 
dom perturbation, and expand it about the unperturbed 
surface R = R s , which is translationally invariant along 
the x-axis [R = (x, r)], so that: 



* (R = Rg) =0, for x < and x > L, 



(3a) 



= -C (R) • for < x < L. (3b) 

Alternatively, the latter boundary condition can be asso- 
ciated to a waveguide surface with a random admittance. 

It can be shown that the matrices of reflection and 
transmission coefficients satisfy the following differential 
equations^: 



dr 



dt 



- (e~ i%L + re fkL ^j v (e~ ikL + e fkL r^j , (4a) 



dL 2 



kL^> { „—ikL i AkL- 



(-' 



(4b) 



with k = diag (k n ) and 



ds 4> m {s)C{L,s)4> n (s) , 



^( S ) = fc- 1 /2 n (r s ). 



dXn (r) 



Or 



it has been assumed that C = C n - The explicit form of 
the differential ds over the cross section (oriented) sur- 
face element ds = rids depends on the geometry under 



consideration. The reflection and transmission intensities 
arc defined by: 



R 



ma |?"mn| i^ran — |^mn| •> 



(5) 



which yield the intensity coupled into the rath outgoing 
channel in reflection and transmission, respectively, for a 
given nth incoming channel. 



B. Field distribution inside the disordered region 

By invoking Green's theorem, the expression for the 
field inside the waveguide (0 < x < L) can be written as: 

*„(x,r) =*° (x,r) 

L 



J dx' I ds'^ n (x',r' s ) 



, i \ dG (x',x;r' s ,r) 

dn' ' [ ' 



where = fc„ 1 ^ 2 x„(r)e lknX . Substituting the Green's 
function 



N 



Go (x,r;xV) = £ (2ik m )- 1 Xm (r) Xm (r')eM— 1 

m— 1 

into Eq. © , we end up with the following expression for 
the scattered field inside: 

^ c (x, !•) = *„ (x,T)-9° n (x,r) 

L 



f2(2ik m )-\ m (r)^^-e ik ^-^ 



x 



7/1 = 1 

>0 _ I. -1 / 2 ,, ik n x 



(8) 



where \P„ = k n x„(r)e nX . Rearranging the inte- 
grand, and handling the phase factor e lk ™\ x ~ x I appro- 
priately by splitting the integral along the waveguide 

L x L 

J = J + J, one obtains: 



a: 



iV 



m=l 

x I jdx'e ik ^-^ jds>* n (x',r' s ) 
{ o 

+ f dx'e~ ik ^ x -^ Ids'V n (x',r' s ) dx ™(? a) \ . (9) 



Factoring out the phase factors defining waves propagat- 
ing right and left, and splitting again the integral of the 
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second term, we get 

Xm(r) f e' fc " 



' A- 1/2 I 277r 1/2 



-ik m x 



2ik. 



1/2 



dn' 



-(.10) 



At this point, we define the local amplitudes of the scat- 
tered waves propagating along the x axis in positive and 
negative directions, respectively, p{x) and t(x): 



p L mn (x) = {2ik 1 J l 2 y 1 

X 

x J dx' I ds'^ n (x',r' s ) 



I „l \ %m( r s) -ik m x' 

On' 



I „r\ ^ Xm( r s ) Ak m x' 

8n> 



T mn (x) — t mn (L) $mn (^^m ) 

x 

x J dx' <j) ds'^ n [x',r' a ) 

so that 

*~(R)= X)C /2 Xn(r) 



(11a) 



(lib) 



x^We^+iWe^ 1 }. (12) 

Then, by differentiating Eqs. 1|11|) . and taking into ac- 
count the boundary condition © in the integrands, with 
the aid of Eq. I|12|l again, a set of coupled differential 
equations for Pmn( x ) an d T mn( x ) is derived: 



dp L 
dx 

d? L 
dx 



e lkx v 



e lkx v 



e- ikx (T+r L )+e ikx p L , (13a) 
e - ik ~ x (T + t l ) + e^ x p L ] . (13b) 



The corresponding boundary conditions satisfied by 
Pmn( x ) an d T mn( x ) & t the end points of the waveguide 
are: 

ptn{ x = Q) = 0,T^ n ( X = 0)=t m n(L)-6 mn ; (14) 
Pmni x = L ) = r,nn(L), T^ n (x = L) = 0. (15) 

C. Numerical calculations 

We have chosen for the numerical simulations the same 
geometry as in Ref. two parallel, perfectly reflect- 
ing planes at z = and z = d with random deviations 
z = £ (x) given by a ID stochastic process with Gaussian 
statistics with zero mean and a Gaussian surface power 
spectrum 



g(Q) = sMaexp (-(Qa) 2 /4) 



(16) 



where S is the RMS height and a is the transverse correla- 
tion length. The corresponding transverse eigenfunctions 
are thus given by 

Xn (z) = (2/d) 1/2 sin (n n z) , «„ = nn/d, (17) 



and the impurity matrix (JjjJ by 



d (&n^m)^ 2 



(18) 



In order to model ID wave propagation in our calcula- 
tions, the waveguide supports only one mode, its thick- 
ness d being such that Q = u)d/(2Trc) « 0.75. Conse- 
quently, all subscripts referring to mode indexes are sup- 
pressed hereafter. 

The linear differential equations for the reflection and 
transmission amplitudes (JIJ are solved numerically by 
means of the Runge-Kutta method; this is done for a 
given realization from L = up to a maximum 

length L = L max (cf. Ref. 0). Then, for a fixed length 
of the disordered segment L, the same standard numeri- 
cal techniques are employed for the system of first-order 
differential Eqs. (|13|) in order to obtain the local reflection 
and transmission amplitudes, with the help of the bound- 
ary conditions l|15(l involving the reflection and trans- 
mission amplitudes [r(L),t(L)], previously obtained. Fi- 
nally, the field intensity is calculated from the incident 
and scattered fields inside [Eqs. (JTJ;) and l(T^)l]: 



I(x)=\*(x,d/2)/V°(x,d/2) 



(19) 



D. Analytical approach: Rapid phase averaging 

To calculate analytically the average intensity, (I(x)), 
inside a one-dimensional disordered system we introduce 
the function 



R(x) 



p(x) 



exp(2i/cx), 



t{x) + 1 

that satisfies the nonlinear equation 

i^ = -2kR(x) + ^[l + R(x)f, 



dx 



(20) 



(21) 



where k is the longitudinal wave number of the propa- 
gating mode (fc = (w/c) 2 — (ir/d) ). The random 

scattering potential, V(x), in the case under considera- 
tion, i.e. in a single-mode waveguide with a randomly 
rough surface, has the form: 



V(x) = ~C(x). 
Then the intensity can be expressed as4 



I(x) = I l-\R(x 



1 + \R(x)\ 2 + 2^(R(x)) 
1-I^)| 2 ' 



(22) 



(23) 
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Obviously, R(L) = r(L) exp(2ikL) where r(L) is the 
total reflection coefficient of a single-mode waveguide de- 
fined by Eq. I|4a|l. ft is convenient, following Ref. 0, to 
introduce two functions, u(x) and ip(x), so that 



R(x) 



lu(x) - 1 
u{x) + 1 



exp [i<£>(a;)] . 



(24) 



Substitution of Eq. (JSJ into Eq. yields 



I(x) 



2In 



u(L) + 1 



u(x) + \/u 2 (x) — lcos(p(x) . (25) 



If the scattering is weak enough, so that l sca t 3> A, the 
random phase, f{x), is uniformly distributed over [0, 2tt]; 
We have verified this assumption through numerical cal- 
culations of the probability density function of <p(x) (not 
shown here), which indeed yield a uniform distribution 
in all cases studied below. Obviously, to get rid of the 
rapid (on the scale of order of A) oscillations of the phase 
one has to integrate (average) Eq. (|25|l over an interval 
£ that satisfies the inequality A <C £ <C I scat- This rapid 
phase averaging (RPA) yields 



I(x) 



2Iqu(x) 
u{L) + l 



(26) 



The two-point probability distribution function, 
P2{ul, L; u x , x), necessary for the ensemble averaging of 
the intensity I(x), Eq. pfi]l. can be also calculated under 
the assumptions that l SC at 3> A and that the scattering 
potential is a ^-correlated Gaussian random process such 
that 



(5V(x)6V(x )} = DS(x - x 



(27) 



Then, the smoothed (RPA) mean intensity distribution 
inside a one-dimensional random system can be presented 
in the form4 

(I(x)) = vrexp[L>(x - i/4)] / dfi f— 

Jo cosh fin 

x exp(-fi 2 DL) (cos2fiDx + &m 2 ^- Dx ^j _ ^8) 

In what follows, comparisons with the numerical calcu- 
lations will be made on the basis of the average (macro- 
scopic) properties, regardless of the (microscopic) details 
of the disorder. Namely, the localization length £, defined 
from (logT) ~ — will be used as matching parame- 
ter, which in this RPA approach is given by £ = D^ 1 . 



III. SINGLE REALIZATIONS: RESONANCES 

First, we identify the roughness and waveguide param- 
eters that lead to the onset of Anderson localization. This 
is done in Fig. [21 by plotting the length dependence of 
(logT) at frequency u> = 0.75 (single mode) for several 



RMS heights 8 and fixed correlation length a/d = 0.2. 
The resulting linear decay is the fingerprint of Ander- 
son localization, the decay rate yielding the localization 
length. The fitted values of £ for each S are included in 
Fig. El 

With the aid of the latter results, we choose a set 
of parameters that ensure the ID Anderson localization 
regime: a/d = 0.2, S/d = 0.05, and L = 1500d « 5.5f. 
We then calculate the frequency dependence (in a nar- 
row frequency range) of the transmission coefficient T(u>) 
for a given realization, as shown in Fig. |3fa). Extremely 
large fluctuations are observed with narrow spikes ap- 
pearing over a fairly negligible background. The latter 
background yields the expected response at typical (high 
probability) frequencies, since (log T(cD = 0.75)) ~ —5.2 
and (T(Q = 0.75)) - 0.062. The low-probability peaks 
in Fig. OUa) correspond to narrow resonances or quasi- 
transparent frequencies at which the transmission coeffi- 
cient can be even 1. 

The transmission in the vicinity of one such transpar- 
ent frequency (fio = 0.75069) is presented in detail in 
Fig-EIb). Note the frequency scale, revealing how nar- 
row the resonance is. By fitting the numerical result to a 
Lorentzian [also shown in Fig.^b)], we obtain the half- 
width at half-maximum T/uio w 2.4 x 10 -6 . Resonances 
behave like high-finesse cavity modes with large associ- 
ated Q-factors (~3x 10 5 ), which may lead to practical 
applications as in random Iasing2*i2iii. 

The field intensities inside the waveguide for fre- 
quencies at the resonance, mid-resonance, and out-of- 
resonance [Q = 0.75069, 0.7506875, and 0.7506, respec- 
tively, in Fig. Olb)] are shown in Fig. where the 
envelope and average of I{x) over rapid oscillations (pe- 
riod ~ n/k) are plotted. The incident mode impinges on 
the disordered segment at x = L propagating from right 
(positive x axis) to left (negative x axis). At resonance 
[see Fig. 0Ja)], high intensity concentration takes place 
over a region around the center of the disordered seg- 
ment of the waveguide (I ~ 200 with a peak of I ~ 400) , 
its particular shape being a characteristic feature of the 
given resonance. The field intensity at the end points 
(not discernible in the figure) is I(x = 0, L) = 1, as ex- 
pected (T = 1,R = 0). At mid-resonance [see Fig. 0Jb)], 
the field intensity distribution maintains its shape, but 
the overall height is decreased by nearly a factor of 2. 
The reflected and transmitted coefficients are retrieved 
at the end points: I(x = L) =| 1 + rexp(ikL) | 2 (enve- 
lope w 2.25 and mean « 1.5) and I{x = 0) = T « 0.25. 

In contrast, an absolutely different behavior has been 
observed away from resonance, i.e. at typical (non- 
transparent) frequencies (or realizations), as seen in 
Fig. 0Jc) . The field energy is not localized, but decays 
from its initial value I(x — L) (envelope « 4 and mean 
w 2) to the exponentially small value I(x — 0) = T ~ 
exp(-i/0. 
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IV. TOTAL, TYPICAL AND RESONANT 
AVERAGE FIELDS 

We now turn to the analysis of the ensemble average of 
the field intensity (I(x)) along the disordered region. Nu- 
merical simulation calculations are carried out for fixed 
u> = 0.75, L, and statistical parameters of the roughness. 
Averages have been done over N = 10 5 realizations, sep- 
arating typical and resonant realizations according to a 
threshold value of the mean transmission coefficient T c . 

FigureEIa) shows (I(x)) for a/d = 0.2, 8/d = 0.05, and 
various values of the disordered segment length L/d — 
1200,1500,2250,3000. (Recall that the incident mode 
impinges on the disordered segment from the right end, 
x = L, which we have shifted to the origin for the sake of 
clarity.) In all cases, the mean intensity decays monoton- 
ically towards the exit of the disordered waveguide, the 
decay rate being smaller the longer is the waveguide (pro- 
vided that L/£ 3> 1). The contribution from resonances 
to the mean intensity, (I(x)) reso , yielding transmission 
coefficients larger than T c = 0.4, is shown in Fig. p>ft>). 
Broad distributions are found with large maximum field 
intensities lying near the center of the disordered waveg- 
uide. The contribution from typical realizations (> 90%), 
(I(x))ty P is plotted in Fig. Efc); a qualitative behavior 
similar to that of the total mean intensity is observed, 
except for a faster decay rate. 

In order to improve our understanding of the physics 
underlying the formation of the field intensity patterns, 
we have replotted (I(x)) by rescaling the x-dependencc 
in units of the disordered segment length L, x = (x — 
L/2)/L. In addition to that, calculations have been done 
for different roughness parameters 8/d = 0.05,0.08 and 
0.1 (fixed a/d — 0.2), by choosing L in such a way that 
the ratio L/£ remains fixed (the corresponding values of 
the localization length £ are given in Fig. |2J). The re- 
sulting (I(x)} are presented in Fig. El The RPA quasi- 
analytical results obtained from Eq. (|28J) are also in- 
cluded. 

Several conclusions can be drawn from the latter re- 
sults. First, (I{x)) exhibits in all cases a universal be- 
havior, depending only on the ratio L/£ regardless of 
the microscopic details of the ID disorder. Actually, 
as shown in the inset in Fig. we have observed that 
universality can be pushed further, so that (I(x)) [with 
x = (x — L/2)/(L£) 1 / 2 ] is a unique function. Sec- 
ond, for moderate and even large the RPA ex- 
pression predicts very accurately the mean field distri- 
bution obtained numerically; a monotonic decay from 
(I(x — 0.5)) = 1 + (R(L)) at the incoming end to 
(I(x = -0.5)) = (T(L)), crossing the value (I) = 1 
through the middle of the disordered segment x = 0, 
and being steeper the larger is L/£. Third, deep into 
the ID Anderson localization regime, £/£ > 11 in Fig.El 
the numerical results reveal a departure from the RPA 
predictions, as evidenced by the shift of the (I) — 1 cross- 
ing towards the incoming end. We have investigated the 
physical origin of this discrepancy by enforcing in the 



numerical calculations some of the assumptions made in 
the RPA approach. First, uncorrelated disorder has been 
used in the numerical calculations, with similar results 
to those for the Gaussian correlation. Rapid phase aver- 
aging has also been carried out at each realization prior 
to ensemble averaging, yielding no significant differences. 
Thus neither finite correlation nor RPA can give rise to 
the observed discrepancy. 

At this point, it is important to emphasize that plot- 
ted in Fig. is the ensemble average of the intensity, 
which is a non-self-averaging (strongly fluctuating) quan- 
tity. To gain insight into the behavior of the field inten- 
sity pattern at different individual realizations, we have 
separated typical and resonant realizations according to 
a threshold value, T Cl of the mean transmission coeffi- 
cient. The contribution from resonances to the mean 
intensity, (I{x)) reso , and (rescaled) (I(x)) reso , yielding 
transmission coefficients larger than T c = 0.4, is shown 
in Figs. EJb) and One can see that the contribution 
from resonances also exhibits universal behavior in the 
form of a broad distribution with a relatively large max- 
ima within the disordered segment, being determined not 
only by the ratio L/£ (as in the case of the total average), 
but also by the the cutoff parameter T c . Actually, from 
the comparison of the curves for (I(x)) reso with differ- 
ent T c in Fig. it follows that T c fixes the position of 
the maximum intensity, whereas the ratio L/£ sets the 
precise value of the maxima. For fixed T c , the maxi- 
mum intensity is higher for larger L/£; namely, stronger 
resonances are needed for longer disorder in order to cou- 
ple the same amount of energy through the system (or 
similarly, to tunnel through a wider barrier). Relaxing 
the definition of resonance (lowering T c ) for fixed ratio 
L/£ (see Fig. [7J|, leads to asymmetrical (I(x)) reso dis- 
tributions with maxima shifted from the center to the 
incoming end of the disorder segment. 

The contribution to the average intensity from typi- 
cal realizations (> 90%), (I(x))t yp , plotted in Fig. |S[c), 
also depends universally on L/£ and T c (not shown here), 
and shows a qualitative behavior similar to (I(x)), with 
a faster decay, as expected. Interestingly, neither (I(x)) 
nor (I(x))t yp decay exponentially, but rather manifest a 
s-like dependence, as mentioned above. This means that, 
no matter how long the realization is [i.e., how small is 
exp(— L/£)], a lengthening of the system (increase of L 
with the localization length kept fixed) leads, surprisingly 
enough, to essential changes in the energy distribution 
inside. Namely, provided that the strength of the disor- 
der is fixed, the longer a randomly disordered sample is, 
the slower is the decay of the intensity [both (I(x)) and 
(I(x))typ\ from the incoming end deep into the sample 
[this is neatly observed in Figs.|Sfa) and Efc)]. In other 
words, the "penetration depth" of both (I) and (I)typ 
into a ID random system is independent of the strength 
of scattering. This effect is, however, dependent on the 
value T c in the definition of (I)typ, as illustrated in Fig. [8] 
with the cutoff decreasing, the slowly decaying part of 
(I(x))ty P near the incoming end (x — L) diminishes, the 
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distribution thus decaying more abruptly. Obviously, the 
longer a realization is, the smaller T c is necessary for the 
transition to take place. Interestingly, the intensity for 
a single typical realization for which T ~ (T) appears 
to decay approximately exponentially / ~ exp(— x/£), as 
seen in Fig. |H1 (its oscillations are smoothed spatially on 
a log scale). This is in accordance with the behavior of 
the average logarithm of the intensity, which fluctuates 
less strongly than the intensity itself and fits very accu- 
rately (log I(x)) ~ — \x — L\/£ (see Fig. 0, revealing its 
self- averaging nature. 

Finally, we have calculated higher-order moments of 
the mean intensity (I n (x)). In Fig. [^1 the numerical re- 
sults are shown in the case n — 2, 4 for some of the dis- 
ordered waveguides considered above. The most remark- 
able feature is the broad, resonant-like shape, revealing 
the increasing (for higher n) influence of (low-probability) 
resonances, with huge field intensities. 

V. CONCLUDING REMARKS 

To summarize, we have developed a formalism to cal- 
culate the field inside surface-disordered waveguides, sim- 
ilar to that of the invariant embedding equations for 
the reflection and transmission coefficients. By applying 
it to 2D single-mode waveguides with planar walls and 
Gaussian-correlated surface roughness, we have investi- 
gated the occurrence of resonances in the ID Anderson 
localization regime, with emphasis on the resulting field 
intensity distribution both for given realizations and en- 
semble averages. 

We have examined the frequency dependence of the 
transmission coefficient T(lo) for different realizations; it 
exhibits well-defined resonance-type behavior inherent to 
the localization regime. This enables us to separate typ- 
ical realizations, characterized by very low (as expected 
from the average (logT) ~ —L/C) values of T and a 
monotonically decaying intensity, from resonances with 
transmission coefficients close to one and extremely high 
intensity maxima (localization) in a region around the 
center of the system. 

Numerical simulation calculations for the mean field 



intensity (I(x)) along the disordered segment of the 
waveguide reveal a universal behavior completely deter- 
mined by the ratio L/£: A smooth decay from the initial 
value of (I) ~ 1 + (R(L)) at the incoming end, to the 
outgoing mean transmitted field intensity (J) ~ (T(L)), 
crossing the value (I) ~ 1 at/near the center of the disor- 
dered segment. For moderately strong disorder L/£ ~ I , 
the quasi-analytical (RPA) prediction H28JI fully agrees 
with the numerical calculations. However, for strong dis- 
order L/£ 3> 1, the numerical results exhibit, unlike the 
RPA result, a shift of the mid-point ((/) ~ I) towards the 
incoming edge. The contribution to (I(x)) from resonant 
realizations (those yielding anomalously large transmis- 
sion above a threshold value T c ) manifests also universal- 
ity characterized by the parameters L/£ and T c : Its shape 
is a broad distribution whose maximum value, which is 
larger for stronger disorder, shifts from the center to- 
wards the incoming edge with decreasing T c . On the 
other hand, we have found that the contribution from 
such low-probability resonances become more dramatic 
in higher-order moments of the total intensity distribu- 
tion. 

The contribution from typical realizations to the total 
average, (I(x))t yp , depends on the cutoff value T c . For 
T c not too small, (I(x))t yp (as well as (I(x))) inside a ID 
random system is slightly dependent on the strength of 
the scattering, and increases with the increase of the total 
length, L, of the system. With decreasing threshold value 
T c , the penetration depth ceases to depend on L and 
{I(x))t yp decays more rapidly. In this regard, evidence 
of the self-averaging nature of \ogI(x)/x is given by the 
behavior of log/(x) for single, typical realizations, and 
also by the result that (log J(x)) ~ —\x — 
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Figures 




FIG. 1: Illustration of the scattering geometry of the surface- 
disordered waveguide. 
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FIG. 2: Average logarithm of the transmission coefficient over 
N = 10 J realizations as a function of the disorder length 
L for surface roughness parameters a/d = 0.2 and 5/d = 
0.0125,0.025,0.04,0.05,0.08,0.1. The localization lengths £ 
resulting from fits to linear decays are shown. 
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FIG. 3: (a) Spectral dependence of the transmission co- 
efficient for a given disorder realization with a/d = 0.2, 
8/d = 0.05, and L/d = 1500 in a narrow frequency range 
showing several resonant frequencies, (b) A single resonance 
is zoomed in and fitted to a Lorentzian (dashed curve, indis- 
tinguishable from the numerical result. 
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FIG. 4: Field intensity along the disorder realization used in 
Fig-C3» (with a/d = 0.2, S/d = 0.05, and L/d = 1500) at (a) 
Q = ud/{2-Kc) =0.75069 (resonance), (b) Q =0.7506875 (mid- 
resonance), and (c) uj =0.7506 (out of resonance, typical). 
To suppress rapid spatial oscillations, the envelopes (upper 
curves) and spatial averages (lower curves) are shown. The 
incident wave is coming from the right end. 




FIG. 5: (a) Spatial distributions of average field intensity 
(N = 10 5 realizations) for a/d = 0.2, 5/d = 0.05, and disorder 
lengths: L/d =550 (curves A), 1500 (curves B), 2250 (curves 
C), and 3000 (curves D). All curves have been shifted to make 
coincide the incoming ends at x = 0. The contributions from 
resonant realizations (with T > T c = 0.1) and the remaining 
typical realizations are shown in (b) and (c), respectively. 
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FIG. 6: Spatial distributions of average field intensity (N = 
10 5 realizations) as a function of renormalized position x = 
(x - L/2)/L for a/d = 0.2 and S/d =0.05 (circles), 0.08 
(squares), and 0.1 (triangles). In each case, several disordered 
lengths are considered according to L/£ = 0.5, 1, 2, 4.4, 8.2, 11. 
Solid curves represent the quasi-analytical, RPA results. In- 
set: the RPA results only as a function of i = (x — 
L/2)/(L0 1/2 - 
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FIG. 7: Same as in Fig. but for the contribution from res- 
onant realizations with T c — 0.4 (hollow symbols, L/£ = 
1, 2, 4.4, 5.5, 8.2) and T c = 0.01 (filled symbols, L/£ = 11). 




FIG. 8: Spatial distributions of the contribution to the av- 
erage field intensity (log scale) from typical realizations with 
a/d = 0.2, 5/d =0.05 (circles), 0.08 (squares), and 0.1 (trian- 
gles), and fixed disordered length L/f = 5.5, for T c = 10 -4 , 
0.004, 0.1, 0.4, and 1 (the latter equivalent to (I(x)), thin 
solid curve). Also included are: exp(log I(x)} (dashed curve) 
and I(x) (spatially averaged in a log scale) for a single, typical 
realization (thick solid curve). 
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FIG. 9: Same as in Fig. © but for the 4th (a) and 2nd (b) 
moments of the field intensity, and L/£ = 2, 4.4, 8.2. 



